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Abstract 

Approximate analytical solution of two dimensional problem for sta- 
tionary Navier- Stokes, continuity and Fourier-Kirchhoff equations describ- 
ing free convective heat transfer from isothermal surface of half infinite 
vertical plate is presented. The problem formulation is based on the typ- 
ical for natural convection assumptions: the fluid noncompressibility and 
Boussinesq approximation. We also assume that orthogonal to the plate 
component of velocity is small. Apart from the basic equations it includes 
boundary conditions: the constant temperature and zero velocity on the 
plate. At the starting point of the flow we fix average temperature and 
vertical component of velocity, as well as basic conservation laws in in- 
tegral form. The solution of the boundary problem is represented as a 
Taylor Series in horizontal variable with coefficients depending on vertical 
variable. 



1 Introduction 

The results of theoretical and experimental study of free convective flows from 
heating objects are widely published and they are very useful to determine 
convective heat losses from apparatus, devices, pipes in industrial or energetic 
installations, electronic equipment, architectonic objects and so on by engineers 
and designers. 

The problem of convective flow development traditionally is based either on 
boundary layer theory or on self-similarity theory [T], [2], [3]. Both methods 
of natural convection heat transfer description use simplifications that allow to 
transform the basic fundamental equations of Navier-Stokes, mass and Fourier- 
Kirchhoff equations. As a main point of the methods is an ordinary differential 



equation introduction which solution relates to a general problem by some con- 
traction procedure. Such description don't give a possibility to investigate some 
details of a flow fleld for example at vicinity of its starting point. It is known 
that the starting point influence is important in many aspects of fluid flow 
development, especially in critical conditions. 

There are some problems of convection for which the boundary layer ap- 
proach fails. One of such cases is an assumption that boundary layer starts on 
a plate, that contradicts experiment. The other is conceptual for the theory: 
the fluid flow is restricted by a conditional boundary which also don't exist in 
reality. Due to this Prandtl introduced integration across a boundary layer to 
determine average values of fluid flow (velocity, temperature) . 

The similarity solution method is based on a speciflc combination of inde- 
pendent variables introduction that allows to get rid of two-dimensionality of 
the general problem. The resulting ordinary differential equation gives an in- 
formation of velocity and temperature profiles that is an advance in comparison 
with boundary layer theory. However such description still is rather qualitative 
than quantitative. For example it doesn't give an answer about flow behavior 
in the vicinity of the starting point in the flow. We would stress that an interest 
to the problem still exists [4], [5], [6]. 

We consider approximate analytical solution of stationary convective fluid 
flow induced by an isothermal vertical surface with the axis y parallel to gravity 
force X is used as a horizontal one. The choice of the coordinate system is 
typical for laminar natural convection simplifications with the following main 
assumptions: fluid is incompressible, Boussinesq approximation, normal to the 
surface component of velocity is neglected [3 [7] . The case of isothermal vertical 
surface is obtained as a limit of horizontal conic with a base angle of the conic 

a = [siini 

Our solutions are build without use neither of boundary layer nor self - 
similarity concepts. We are looking for the solution as Taylor series in x and 
derive a system of equations for coefflcient of the expansion. Each cutting of 
the series results in polynomial approximation in x and flnd a closed system 
of equations in y. We restrict ourselves by upper half plane [y ^ 0). Such 
approach needs a formulation of boundary conditions in the starting point of 
the flow. We take into account the physical character of the flow: namely its 
natural convection origin that means average velocity is zero at j/ = level. For 
the formulation of boundary conditions we base on conservation laws in integral 
form. 

2 The basic equations 

Let us consider a two dimensional stationary flow of incompressible fluid in the 
gravity field. The flow is generated by a convective heat transfer from solid plate 
to the fluid. The plate is isothermal and vertical. In the Cartesian coordinates 



2 



X, y the Navier-Stokes (NS) system of equations have the form [T]: 



OX oy J oy \ oy^ ox 

W^ + W^\=-^+pv(^^^ + ^^^\ (2) 
^ 9x ^ dy ) dx \ dy'^ dx^ J 

In the above equations the pressure terms are divided in two parts p — po + p. 
The first of them is the hydrostatic one that is equal to mass force —gpoo, where 
P = Poo {1 — b{T — Too)) = Poop' is the density of a hquid at the nondisturbed 
area where the temperature is Too • The second one is the extra pressure denoted 
by — Vp.The part of gravity force gb{T — Too) arises from dependence of the 
extra density on temperature, 6 is a coefficient of thermal expansion of the 
fluid. In the case of gases b = \ 'Ij^] = T^^.The last terms of the above 

equations represents the friction forces with the kinematic coefficient of viscosity 
v. 

The mass continuity equation in the conditions of natural convection of 
incompressible fluid in the steady state has the form:. 



dW^ dWy 
dx dy 



0. (3) 



The temperature dynamics is described by the stationary Fourier-Kirchhoff 
(FK) equation: 

dx ^ dy \ dy"^ dx"^ J 

where Wx and Wy are the components of the fluid velocity W that are shown 
on the Fig.l; T, p - temperature and pressure disturbances correspondingly and 
a is the thermal diffusivity. From the point of clarity of further transformations 
we use the same scale I along both variables x and y. We will return to the 
eventual difference between characteristic scales in different directions while the 
solution analysis to be provided. After introducing nontraditional variables: 

x' = x/l,y' = y/l,T' = (T - T^)/{T^ - Too), 

P' = P/POO, = Wx/Wo, Wy = WyjWo (5) 

we obtain in Boussinesq approximation (in all terms besides of buoyancy one 
we put p~ Poo)- 



, dW' , dW' 

W + W 

dx' y dy' 



L (T' + 1) 



Pc 



dp' 



PooW^ dy' 



W, 



,9W'x , ,,„dW'x 



dx' 



W, 



Po 



dy' 



Tiv^d^ 



'd^W' d^w'S 



Qyl2 



dx'^ 



d'^WL d'^WL' 



dy" 



dx'^ 



(6) 
(7) 
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and FK equation is written as 



„„9T' ^.or ,(d^T' d^r 

W h W = a 1 

"9a;' y dy' \ dy'^ dx'^ 



(8) 



-^jy2 — i, Hs a characteristic linear dimension 

and Wq is characteristic velocity. If 



r 



then a' = Pt , ly' = 1 and L = Gr, where 

ijrr ' — " ^ ■ 



(9) 



(10) 



(Grashof number) such link will be used in final solution. After cross differen- 
tiation of equations ^ and ([t]) we have: 



,dW' ,dW' , , 



dx' 
d_ 



dy'' 



dx'' 



" dx' y dy' \ dy'^ 



dx''- 



(11) 
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The FK equation rescales as 



^,ar: arx /a^ 

^ dx' y dy' ) \ 9a;'2 / 

Next we would formulate the problem of free convection over the heated vertical 
isothermal plate x = 0, y G [0, Z), dropping the primes. 



3 Method of solution and approximations 

Next aim of this paper is the theory application to the standard example of a 
finite vertical plate. In this case we assume the angle between the plate and a 
stream line is small that means a possibility to neglect the horizontal component 
of velocity of fluid, denoting the vertical component as W (y, x). It means that 
the number of equations: (11), ( [l2| ), (|3| exceeds the number of variables: T 



and W , hence one of the equations should be excluded. Having in mind mass 
conservation account in integral form, we exclude the equation ([5]). The mass 
conservation law introduces boundary conditions as a link between mass law 
flow on levels y — Q and y ~ L (see the section on boundary conditions). Con- 
sider the power series expansions of the velocity and temperature in Cartesian 
coordinates: 

W {x, y)=j{y)x + a (y) x^ + p (y) x^ + >f(y)x'*..., (13) 

T (x, y)^C{y)x + A (y) x^ + B (y) x^ + F (y) x\... (14) 

According to standard boundary conditions on the plate we assume that the 
both functions tend to zero when x ^ 0, we choose for a calculation convenient 
the zero value for nontraditional temperature ([5|.It means that the value of 
T{x,y) outside of the convective flow tends to —1. 

We would like to restrict ourselves by the fourth order approximation for 
both variables that means we neglect higher order terms starting from fifth one 
(see the Fig.l where the area of the approximations validity is marked as dashed 
one). 

As it will be clear from further analysis we should consider the functions: 
<^{v) T P [y] 1 C (y) and B (y) as variables of the first order, while 7 (y) and F (y) 
to be the second one. From the relations that appear after substitution of ([T3| 



and (141 into m and (??) it follows that A{y) = 0. 



Substituting expressions (13) and ( [14| into the equations (??) and (??) 
taking into account the assumption that Wx = and Wy = W yields 



d_ 
dx 



W— -G (T+1)- f — + — 



= 0, (15) 



From ([l6| one found that A{y) — and F (y) = 0. Finally from both equations 
(15), (16) we obtain the system of equations for the coefficients B{y), C {y), 
a [y) and /3 {y) : 



&B [y] 



d'C{y) 



= 0, 



dydy 
dC{y) d'B{y) 



i{y) 



dy dydy 
6/3 {y)-GrC{y)=Q, 
dj (y) d'^a (y) 



dy 



dydy 



0. 



(17) 

(18) 
(19) 
(20) 



The first two (17), (18) arise from FK equation and the rest of them are from 
the NS one. The system of equations is closed if 7 {y) = const = 7. It means 
that the number of equations and the number of unknown functions is the same. 
Finally, in the first approximation the velocity and temperature are expressed 
as: 

W{y,x)=^x + aiy)x^+(3iy)x\ T (y, x) ^ C (y) x + B (y) x\ (21) 



From (1201) one has 



a{y) = Ciy + C2. 

6 dydy ' 



(22) 



From (FtI) it follows that B (y) = -i^^^, hence (|18|) goes to: 



1 d^Cjy) _^ ^ dCjy) ^ ^ 

6 dydydydy dy 



The equation ( 19 ) reads 



Hy)^~C{y). (24a) 

The equation (|23]) is the ordinary differential equation of the fourth order, there- 
fore its solution needs four constants of integration, denoted as: C(0), C/(0), 
C"(0) and C""(0). These constants depend on two parameters Ci and C2, which 



enter the coefficients of the Eq 23 The function C{y) defines the rest functions 
/3 (y) and B (y) via above relations. It means that we have six constants deter- 
mining the solution of problem and we need also six corresponding boundary 
conditions at j/ = 0. 
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Sketch of the natural convection flow for air (Pr =0.7). 



4 Boundary conditions formulation 



At the starting horizontal edge of the vertical plate the bouncy force yet does 

not act therefore we assume that the mean vertical component of the velocity 

is negligibly small. Thus we neglect a fine velocity structure of transition area. 

Similar the temperature values are averaged in the vicinity of the boundary edge 

point and taken as value -1 (temperature of incoming from the bottom flow). 

Then let us define the mean value of a function F{x) (velocity or temperature in 

1 

the interval x e [0, 1] (Fig. 2) in the nontraditional form as: F = J F{x)dx. In 



dimensional form the interval of averaging length is I which hence is defined by 
this integration boundary. Let us remind that scale I is connected with special 
(local, horizontal) Grashof number (10). In natural convection problem 



theory vertical velocity component and temperature values in surrounding of the 
heated plate are: Wy(0,0) = 0, T(O.O) = —1 according to the above notations. 
Taking it into account we put: 



W{Q,x) = 



{^x + a (y) a;^ + /3 (y) x^^ dx 



27+3«(0) 



1 



/3(0)=0 
(25) 



7 



therefore due to ( 22 1 we arrive at the first boundary condition: 

3 



a (0) = C2 



r/3(0) 



,7 



L 



cm 



,7 



(26) 



Nontraditional temperature of the fluid at the lower half plane, according to 
above, is —1. In the analogy to the condition for velocity we a ssu me that average 
temperature at the limit y = is also —1, therefore from (21) T {y — 0,x) = 



(C {Q)x + B (0) x3) dx = IB (0) + (0) 



relation arises: 



Plugging B (0) 



B(0) + 2C(0) 



-1. Then the next boundary 



(27) 



dydy 

second boundary condition: 

1 
6 



1 r^ii 



C"(0) from (??) into (271 gives the 



C"(0)-2C (0) 



(28) 



The third and forth boundary condition arize on the level y = L, where the 
fluid lose the contact with the heated plate. At this point we suppose that 



and 



dT. 
dy 

dW 
dy 



y=L 



= 



\y=L 



= 



(29) 



(30) 



consequence of the equation ( 29 1 is 



Both conditions ( 29 ) and ( 30 1 have a transparent physical meaning: the tem- 
perature and velocity does not grow from the level y = 0. The immediate 

(31) 
(32) 



dC{y) 
dy 



\y=L = 



and 



dB{y) l d^C{y) 



dy 



6 dydydy 



From ( 30 1 one can see that 



da (y) 



y=L — Cl — 
_ dCjy) I 



(33) 



dy 

and, according to (24a), ^g^^^ \y=L = ~^ °'"a'y' \y=L = coincide with (31) 



Looking for the rest boundary conditions let us apply conservation laws of mass, 
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momentum and energy for derivation of boundary conditions at y = 0. The first 
of them is the conservation of mass in steady state in two dimensions looks as : 



/ 



pV^ ■ ltdS = (34) 



where S is a closed surface of a control volume, which cross-section is shown in 
Fig.2. 

The second and third of them are produced by the law of momentum vector 
conservation, which is represented by Navier-Stokes equations ^ and ([t]). Let 

us denote |f = p^^gb (T - T^)-pW^^-pWy^+piy + ^) = Ey 

and If = -pW,^ - pWy^+pi^ + ^) = E,, then one have the 

equation: 

dEy dE^ 

that we use as a basic one. The integral version of this equation is obtained by 
Stokes theorem: 

rotn~tdS = / ~td~t = 0, (36) 



where s is the boundary curve of the integration area S (Fig. 3). 

One of them has the form ( 36 1 . If we specify the curve of integration s as the 
border of the control area S and denote the tangent projection as Et = ■~t 

we have 

' Etds = 0. (37) 



The next one is the first equation of Navier- Stokes ^ : 

The excess pressure is the case of natural heat convection problem to be con- 
sidered is small and proportional to the excess of the temperature: 



dT 



{T-T^)^ p^B! {T-T^). (39) 



In the case of the gas one can use the ideal gas equation of state in dimen- 
sional form p = pR' [Too + {T - T^)] = Po + P ~ PooR' [Too + (T - T^o)] where 
Po — pooR'Too and W = R/p. {p is a molar mass) is gas constant hence 
P = PooR' {T — Too) ■ Going to the nontraditional form, plugging (39) to (381 
yields 
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where we introduced the new nontraditional parameter 



H 



{ghv') 



3 



(41) 



where T^j — Too = W^o and Gr are defined by ([9| and (10) . The viscosity 
coefficient v' here is dimensional. 

The next boundary condition is connected with the conservation of energy in 
a control volume V (area S with unit width see Fig. 3) arises from FK equation 



( 12 ) by integration over the volume. 



Pr 



V 



V 



dV 



gradr)7td5 (42) 



Let us now formulate the rest four boundary conditions on the base of the 
conservation laws. 



According to mass conservation law (34 1 the mass flux don't depend on y. If 

1 



we recall that the mean density is approximately constant, then 



,Wdx 



0, or 



^ {-fx + a {y) x^ +13 {y) x^) dx 



= 0. Integration and differentia- 



tion yield the third boundary condition: 

-lda{y) , ldp(y) 



3 dy 4 dy 



(43) 



y^O 



Plugging the explicit expressions for a (y) ( 22 ) and /3 (y) ( 24a I into ( 43 ) yields 
It gives a link between the constant of integration Ci and the boundary value 



of = c'{y) at the point of j/ = 



dy 



(44) 



Let us apply the conservations laws in integral form to the control volume 
V (see Fig. 3) based on the interval y e [0,£] and x G [0,1]. According to 
our main assumption about two-dimensionality of the stream we neglect a 
dependence of variables on z coordinate. Going to the mom entu m conserva- 
tion let us substitute the components of the vector : into (151 or plugging 

Gr (T-f l)-(i 



dp 
dy 



dy 



dy"^ 



into (37) results in (see Fig. 2): 



)Etds = / Eyl^^ody = 0, 



(45) 
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Fig.3 Notation of coiisei-\'atioiis laws in integral fonii for 
control domain of the iM oblem. 



because the rest parts of the integral are zero according to assumption that 
velocity component = is zero Ej; = 0. 

The left side of the equation ( 45 1 according to ([5| and the fact that the 

velocity W vanishes on the vertical plate surface has the form : J Ey\x=ody = 



e e 



GrE + 



therefore Gr 



^^^^l^U^ody = Gre + 2ja{v)dy = 0, 





lim.^o^ /a(y)d2/ = -2a(0). 



After substitutions of a{y) ([22]) and C2 from (26) one have two relations. 
First one links the integration constant C2 with parameters of the problem 



C2 



while the second one gives the boundary value of the basic function C (0) 

C(0) = ^7-4. 



Plugging the result into (28) gives: 

C"(0) 



144 



■7 - 24. 



(46) 



(47) 



(48) 



The momentum conservation in differential form is represented by (40) that 
we use as a direct source for boundary condition formulation. To determine 

the fifth boundary condition let us evaluate the mean value of the momen- 
tum component balance (shown in above equations) that yields the equation: 
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f {^x + a (y) x^+(3 (y) x^) [^x^ + '-^x^) ^ 
Gr {{C {y)x + B (y) x^) + l) + H [^x^ + 



( 



dydy 



dydy 



which gives the following at the point y = 



x^^ ~i6l3{y)x + 2aiy)) 



dx — 0, 



^ Gl + ^Grj]C'iO) + il2^-2Gr)^0. (49) 



The last (sixth) boundary condition we found from the energy conservation 



equation (42) for the control volume V shown in the Fig. 3 
V s 



(50) 



Taking into account that integral with respect to z gives the factor 1, yields 



Pr 



After substitutions W and T ( pT| ), the left side of ( |50| gives 

VtO 504/ dydydy ' ^ ^ V 120 



Pr (IS; ~ 15) ^'^dy^ , while the right side after transi- 



tion e ^ tends to: 



1 d'Bjy) 1 d'Cjy) 
4 dydy 2 dydy 



y=o 



C(0), 



Plugging the expression for B {17) and next using the equation (23 1 for the 
basic function C{y) gives: 



1 d'C{y) _ 1 d^Cjy) 

2 dydy 24 dydydydy 



C{0) 



G._Pic'(0) + f7- 
The complete equation ( 50 1 gives the last boundary relation: 

7 Gr 



' 70 - 504 j + li) - 1^ ) ^ -^^'(0) + 



(51) 



5 The explicit form of boundary conditions for 

C(y) 



Let us solve the conditions (49 1 and (51 1 with respect to the first and third 
derivatives of C (y) at y = 0. 



The relation ( 49 1 yields at y = : - ^ HC" {0) + {^H - j^Gj + ^Gr^/) C (0)- 
127 - 2Gr = 0. 

" jGrj) C'iO) + ^7 - f a into (|51 1 

' dCjy) ' 
dy 



Plugging C""(0) = (12 - ^Gl , ,40^. 



gives the equation for 



y=o 



12 



V70 



504 



)[( 



12 - 



70H 



140H 



C'(0) 



2887 



48G^ 



' V 120 



M- 



Gr M 
8 



= 



6O7 I 8_ 
PrG^ Pr ^ 

Solving this linear equation yields: 



M = 



Pr G^ 



15 



504 



Gr) 



1 (^2 1 
1680^1- ' 



13 
3360 



Grl) 



(52) 



where H is defined by (41 1 The result allows to evaluate the third derivative at 
y — 0. Let us list the boundary conditions for C{y) : 



C" (0) = f 12 - 



,C' (0)-M,C" (0)-i^7 



70i? 



13 



lAQH 



^ 1 288 



2A, 
48 



(53) 



6 The solution of the boundary problem for C(y) 



Let us rewrite the equation (|23|) substituting the expressions for Ci (|44|) and C2 

G^Pr 

dC{y) 



(61) and introducing the Rayleigh number Ra = Pr . We have 
d^C{y) 



dydydydy 



3i?a 



M 



-y-1 



dy 



= 0. 



(54) 



As it is seen from estimation of M (Gr) one can neglect the term ^y << 1 
in comparison with unit within the interval y € [0, 10] if the Grashof number 
Gr > 500 . Such approximation yields the equation with constant coefficients 
with solution as linear combination of exponents: 



G{y) =^Aexp[kiy], 



i=0 

where ki are roots of the equation 

- 3Rak = 0. 

It means that 



ko = 0, ki 



'iRa = S,k2 



, «'2,3 



( 1 + -^^ 



Plugging (57) into (58), taking into account reality of G{y), results in: 

s \/3 %/3 

Giy) = Ao + Ai exp[sy] + exp[ — ^](Bi cos[— sy] + B2 sin[— sy]). 



(55) 



(56) 



(57) 



(58) 
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The second term Ai exp[sj/] in the above expression exponentially grows for 
positive s, therefore the value of Ai should be analyze in details especially for 
large values of s. 



Boundary conditions ( 53 ) subsequently give 



C(0) = Ao + Ai + Si = ( ^7 - 



(59) 



C"(0) = sMi-Bi; 



Bo 



= M, 



(60) 



(61) 



C""(0) = s^(Ai+Si) =12- 



— Gl + —Gri\M+—^--Gr, (62) 



Solving the system of above equations with respect to Ai and Bi yields 



B, 



1 



M 24 



Bi 



Ai = 



1 M 



1 M 
3 s 

-4^- 



192 
96 



32 
16 



144 



Gr 



48 7 

Gj- 

48 7 



1 M 
1 M 



G„ 



210H s3 



An = - 



(l2-^G2 + ^G.7)A/+^7- 



H 



G ^ 



12^ 

Gr 



13 M 



mH s- 



(63) 

iGr 

(64) 

(65) 
(66) 



The parameter 7 in all above expressions is still not defined. The physical mean- 
ing of it relates to the velocity profile, namely it is the angle of inclination on 
the plate pTj ) . In the literature the profiles are given either directly from exper- 
imental data [H [6] or from theoretical models f2^ . For example in the book [2] 
the velocity dimensionless counterpart f7 = ^ (in our paper is denoted as W) 
in appropriate units is plotted against the function 77 for different Prandtl num- 
bers. The velocity U = '^ 2v^TJr" function of y which is normal directed 
to the plate and x is along the plate (opposite to ours notations). Numerically 
differentiating the expression for U with respect to y on the level x = I and 



and the derivative ^|j7=o 



taking into account the expression (21) we derive the link between constant 7 

C: 



7 



^Gl'\. 



(67) 



The alternative approach to the solution of the equation ( 54 ) is the power series 



expansion. To investigate the parameter values case of order one we use the 



Frobenius method. According to the method let us solve the equation ( 54 ) as a 
power series of y: 
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C(y) = ao + aiy + a2y^ + 03?;^ + 04?/* + a^y^ + ... (68) 



The equations (54) and (68) yield the Unks: 04 = ^aiRa, 05 — ^a^Ra 



j^MaiRa, ag = -^a^Ra — -^Ma2Ra- The boundary conditions give: 

ao = C(0) = ^7-4, (69) 

-i = ^l.^o = M, (70) 

I d^Cjy) 72 

1 a^C (y) , / 1 ^2 13 ^ \ , , 48 8 ^ 

a4 ^Mi?. aa = 4(^7- 12) i?.-4M^i?., 
^« = To " 1680i7 - + 3360i7^ j + T " ^ 

7 The example of air 
7.1 A solution of first type 

For representation of resuhs we choose air as a fluid with typical parameters Tav 
= 20°C,p = 760 Tr, Pr = 0.7 and assume $ = lOK. The physical table values 
in such conditions are following: a = 2. 142 x 10^^ — ,1/ 1. 506 x 10"^ — , i? = 



287.06 = 3. 43 X lO^^^^.For such values the parameter H (Al) is: 



i? = 9.606 4 X lO^G^ (72) 

In the paper [5]. an expression of vertical velocity u{x,y) is derived by the 
method of similarity solution. For the example of Prandtl Number Pr = 0.7 



(air) the coefficient is approximately C = 3/5. Hence the parameter 7 (67) is 
expressed as: 

^ = ^G?/4 (73) 



Our solution depends on the parameter M via equation (54) coefficient. The 



condition in which the coefficient is approximately constant therefore needs the 



parameter M estimation. For the estimation let us express the M (52) as 
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a function of the only parameter Gr with the rest physical parameter values 
account. 



M = - 



1.35 X 10"^/G;- 37634.Gi? + 1.045 x IQ^G} 



3360G^ 



8.59 X W 



1. 05 X lO^G,. + 1. 3 X lO^G? - 11. 6GI + 33. TG^^^ + G/" 



(74) 



M 

2.0 


















1.5 


















1.0 


















0.5 


















0.0 


100 


200 


300 


400 


500 


600 


TOO 





Fig. 4 The dependence of the parameter M against Grashof number Gr 

On the base of the result (see Fig. 4) we choose the Grashof number range 
M/4 < 0.1 for example evaluations of velocity and temperature profiles as a 
function of x and y.li means that approximately G^ > 200. 

For the given values of local Grashof number equal to Gr — 200, 500 and 



1000 the evaluation of the parameter 7 by the expression (67) gives: 7 = 45 



127, 89. 721 and 150. 89. Next the calculation of the solution parameters yields 
the table: 

s = ^i-Gr- 0.7 = 7. 488 9(200), = 10. 164(500), = 12. 806(1000); 

M = 0.297 05(200), 6.1589 x 10-^(500), 1. 2004 x 10-^(1000); 

^0 = -1.3009(200), -1.8475(500), -2. 1894(1000); 

Ai = 6. 652 3 X 10-2(200), 8. 191 1 x 10"3(500), -4. 281 5 x 10-^(1000); 

Bi = -5.803 6 X 10-2(200), -7.487 3 x 10-^(500), 4.3501 x 10-^(1000); 

B2 = -6.452 2 X 10-2(200), -6.786 x 10-^(500), 8. 536 5 x 10-^(1000). 



Analysis of the dependence of the coefficient Ai on Grashof number shows 
its change of sign (see the values for Ai for Gr = 500 and 1000). Plugging the 
expression for M {74) into the formula for Ai (65) and its numerical evaluation 



allow to plot the dependence ^i(Gr).The corespondent plot for the case of air 
is given below (Fig. 5) and allow to determine a critical value of the number 
Gr = Gr cr for which the characteristic linear dimension of the problem has the 

concrete value / — ( ^^^"^^ 
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Fig. 5 The dependence of the Ai coefBcicnt as the function of Gr in the vicinity 

of the critical value G^^cr- 

For Ai w 0, form the plot we estimate Gr = 714. 275, and next calculations give 
7 = 150. 89, s ^ 11. 447, M ^ 2. 916 2 x lO^^, 

Aq = -2. 030 7, B2 = 3. 074 8 x 10"^ Bi = 2. 333 2 x 10"^. (75) 
Plugging the all the values for four cases {Gr = 200, 500, 714.275 and 1000) 



into the solution (581 we plot the resulting curves (Fig. 6). 




Fig. 6 The dependence of the basic function C on y for Gr=200 (red), =500 
(black), 714.275 (blue) and =1000 (green). 

As we announced in the introduction this study is focused on the problem of the 
starting point of the stationary flow in the vicinity of the point x = 0, y = 0. 

For further illustration of the velocity and temperature profiles we choose 
the critical value of Grashof number Gr.cr = 714.275. Pugging the parameters 



values (75) into formulas for the velocity and temperature (21) we obtain the 
profiles of velocity (0, x) = jx+(-^^^ + (-^^ (g^^^^)) x^and 

W{l,x) = 7a; + (%^Mi;- x^ + (-%^C(1)) x^ on the Fig.7 and 

temperature obtained in similar way on Fig. 8. 
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Fig. 7 Profiles of velocity along x on levels y ~ (black) y=l (red) 
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Fig. 8 Dependence of temperature on x at the level y=0 (black) and y=l (red). 
7.2 The solution of the second type for C(y) 



Let us recall the basic relations that define a solution of Navier-Stokes and 
Fourier-Kirchhoff equations via the power series in y. Plugging expressions for 
coefficients 7, a (y) , /3 (y) , C(y) and B{y) in velocity and temperature expan- 



sions (21 ) one have 



W {y, x) = xj + . 



> Gr 



1 



-My - 1 



^'^C{y) 



T{y,x) = C{y)x- 



ia^£(y)^3 
6 dydy 



(76) 



(77) 



where: C{y) = gq + aiy + a2y^ + a^y^ + a4y^ + a^y^ with the coefficients from 



(??), ( [70| and ( [71) ). According to the physical essence of temperature profile 
the function C{y) should be negative. It means that at least the coefficient 
C(l) =00 + 01+02 + 03 + 04 + 05 < O.In such condition of < y < 1 the 
function C(y) < 0. Plugging the expressions for 7, H and M to C(l) yields the 
complicated function of Gr- The result is illustrated in Fig. 9. 
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Fig. 9 The enlarged fragment of C(l) the vicinity of upper critical Grashoff 

number. 

From the plot Fig. 9b it follows that the values of C(l) are negative if Gr > 
Grfl ■ For an illustration we investigate the approximate solution with the local 
Grashof number value equal to Gr — 676.5, for which C(l) = —0.455 29 that 
lays inside the interval (0, —1) which guarantee the negative temperature values 
for all X € (0,1). For the plotting we need the following general parameters 
([74| for Gr = 676.5. The values are: 7 = 112. 56, = 7.403 x 10^ M = 
3.2974 X lO^^and the coefficients ( |70|7l[ ) of the polynomial C{y): oq = -2. 



3.3 X 10" 



, 0-2 



0.0207, 03 = 6. 594 8 x IQ-^, 04 = ^MRa 



003 4, ai = M 

1.9519, 05 = -0.482 02. Substituting these values into Q gives: C{y) 

^y3_0.182 32y2 + 2.9162xl0-2y^ 



560 5y^ + 1. 822 6/ + 5. 832 5 x lO^^y^ _ 0.182 32y2 + 2. 916 2 x 10"^- 2. 030 4.1n 



similar way we derive the expression of C(y) for the critical GraschofF number of 
the previous solution Gr,cr = 714.275. Below wc plot both curves on the Fig. 10. 



0.0 0.1 0.2 03 0.4 oj 




Fig. 10 The power series representation for C(y) at GrashofF number Gr= 676,5 
(red) and Gr,cr = 714.275 (black). 

The reason of discrepancy between the curves is explained mainly by the de- 
pendence of the coefficient 02 on Graschoff number (Fig. 11). 
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Fig. 11 Dependance of a2 on Grashoff number. 

As one can see the value of 02 change the sign at about Gr = 672 

Finally let us substitute the function of C(y) for the Gr,cr ~ 714.275 into the 
expressions for the velocity and temperature and plot them at the levels y = 
and y ^ 0.5(see Fig.l2 and Fig.13). 




Fig. 12 Velocity profiles for the second solution at the y=0 level (black) and 

y=0.5 (red). 
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Fig. 13 Temperature profiles at levels y=0 (black) and y=0.5 (red). 



8 Analysis of the solution and conclusion 

According to the order of power series expansion in x the expressions for velocity 



and temperature in the form of ( 76 1 , [71] satisfy the N-S and the F-K equations 
up to the term . Hence the solution is valid in a stripe adjacent to the surface. 



The equation for the coefficient function C ( 54 ) may be applied on the interval 
y G [0, oo) but the approximate solutions we study here give satisfactory results 
in a vicinity of starting point of the flow y — Q. 

The present theory allows to include higher terms of the expansion hence to 
obtain results that are valid in more wide range of coordinates. Our calculatins 
were performed for the air case that corresponds to the choice of state equation 
in ideal gas form and Prandtl number Pr — 0.7. There is an obvious possiblility 
to extend our results for other fluids, for example as in [S], [5]. 
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